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Abstract 

A recent experimental observation has shown efficient transport of Mega Ampere of electron 
currents through aligned carbon nanotube arrays [Phys. Rev Letts. 108 , 235005 (2012)]. The 
result was subsequently interpreted on the basis of suppression of the filamentation instability in 
an inhomogeneous plasma [Phys. Plasmas 21 , 012108 (2014)]. This inhomogeneity forms as a 
result of the ionization of the carbon nanotubes. In the present work a full 1D3V Particle-in-Cell 
(PIC) simulations have been carried out for the propagation of relativistic electron beams (REB) 
through an inhomogeneous background plasma. The suppression of the filamentation instability, 
responsible for beam divergence, is shown. The simulation also confirms that in the nonlinear 
regime too the REB propagation is better when it propagates through a plasma whose density is 
inhomogeneous transverse to the beam. The role of inhomogeneity scale length, its amplitude and 
the transverse beam temperature etc., in the suppression of the instability is studied in detail. 
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I. INTRODUCTION 


The generation of relativistic electron beams ^EB) through the interaction of a high 


power laser (I> lO^^W/ crri^) with a solid target jl|, and it’s collimated propagation through 
a plasma is a topic of extensive studies due to its applications in cutting edge technologies, 
su* a. fast ignitjon acheme (FIS) S of laaaa fua.o„. — particle accelecafion Q. faat 
plasma switches [5| and laser induced radiation sources [6| . 

The transportation of an electron beam which carries a current much higher than the 
Alfven current limit / = (mc^/e )7 = IT'jbk A, (m is the electron mass, e is the electronic 
charge, c is speed of light and 7 ft is the Lorentz factor of the beam) is not permitted due to 
its own strong magnetic field. However, in plasmas a spatially overlapping return shielding 
current by the background electrons inhibits the magnetic field generation and permits the 
propagation of beams with such high forward currents. The combination of the forward 
and return currents in plasmas are, however, susceptible to several instabilities viz., two- 


stream B, Weibel B, filamentation joj] and oblique modes [1^. The two-stream instability 
is a purely longitudinal electrostatic mode which generate density stripes perpendicular to 
the beam flow direction. The Weibel instability is purely transverse electromagnetic mode 
which occurs in the system due to the temperature anisotropy. This instability generates 
a magnetic field in the unmagnetized plasma system. A similar instability which generates 
a magnetic field in unmagnetized beam-plasma systems is known as the filamentation in¬ 
stability or beam-Weibel instability ju]. These instabilities have a detrimental influence 
on the propagation of the relativistic energy electron beam through a plasma and hence 
put practical constraints. In order to suitably utilize the relativistic electron beams these 
constraints on the propagation of beam through a plasma needs to be overcome. A lot of 
experimental, theoretical and simulation work has been devoted to control and suppress 
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these instabilities |l2l-ll f)l| . The effect on the filamentation instability of beam density m, 


n 
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beam velocity (l^, transverse temperature (l7|] and collisions m has been widely studied. 


A recent experimental [19!] result, has shown many fold improvement in the propagation of 
a hot relativistic electron beam through an array of carbon nanotubes (CNTs). An expla¬ 


nation of this has been put forth by Mishra et al. [20fl suggesting that the inhomogeneous 
plasma created by the ionization of the CNTs by the front of the beam is responsible for 
the stabilization of transverse instabilities, thereby aiding the collimated propagation of the 
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beam through longer distances compared to a homogeneons plasma. Here we carry out 
the PIC simulations to stndy the role of inhomogeneity with arbitrary amplitnde, on the 
propagation of electron beam through plasma. 

The manuscript has been organized as follows. Section II contains the details of the model 
configuration and governing equations. In section III, the analytical results are presented. 
The details of PIC simnlations are given in section IV. The observation of results from PIC 
simulation are presented in section V and it’s interpretation are given in sec.VI. In last sec. 
VII, we conclude our results. 


II. MODEL CONFIGURATION AND GOVERNING EQUATIONS 


The model configuration chosen for our studies has been shown in Fig. [H The beam and 
the return current are chosen to flow along the ±.y respectively. The ion density is chosen 
to have sinnsoidal inhomogeneity along x riding on a constant density of uq. The amplitude 
of the inhomogeneity is e and the spatial variations are characterized by a wavenumber 
kg = 27r/As as shown in Fig. [H The ions are considered to be at rest at such fast electron 
time scales and provide for charge neutralization as a backgronnd species. The density, 
velocity and temperature are denoted by v^, and with suffix a standing for b (beam) 
and p ( background) plasma electrons respectively. The combination of the beam - plasma 
system can be represented by a conpled fluid Maxwell equations, with electrons contribnting 
to forward and retnrn currents treated as two distinct fluids. The normalized governing 
equation in such a case is: 

dua 


dt 

dPa 

dt 


+ V • (riaV^a) = 0 


+ • Vpa ^-(E + 'iTc.xB 


VP. 


rin 


dE -- 

— = V X 5 - EaJa 
dt 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 


with Va = Pa /(1 + Pa'^)K Xx = —'f^a'Pa- The pressure is provided by the equation of state. 
In the above equations, velocity is normalized by speed of light c, density by no, frequency 
by ujq = /rrif, and electric and magnetic field by Eq = Bq = meCUo/e where rrie is 

electron rest mass and e is electron charge. 
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In the equilibrium there is no electric and magnetic held, so there is complete charge as 
well as current neutralization. This is achieved by balancing the forward and return electron 
currents at each spatial location. The total charge density due to electrons is balanced by 
the background ion charge density. The background plasma has been chosen cold (Top = 0 ) 
in all our analytical as well as in simulation studies. The prohle of transverse temperature 
Ta in beam is chosen in such a way that gradient of pressure is zero in equilibrium. The 
temperature parallel to beam propagation direction is chosen to be zero. Thus, for an 
inhomogeneous beam plasma system considered by us in this work we have the following 
conditions for equilibrium: 

nQi{x) = nob{x) + nop{x) ( 5 ) 

S„no„(a;)noa = 0 (6) 

In equilibrium beam pressure To6 is chosen to be independent of x. This is achieved by 
choosing the beam temperature To&(x) to satisfy the following condition 

Pob = Tob{x)nob{x) = constant = k. ( 7 ) 

Tob{x) = k/nob{x) (8) 

The suffix 0 indicates the equilibrium helds. This choice is very artificial construct. However 

it has been chosen to satisfy the equilibrium conditions so that comparisons with simulation 
can be made. We linearize the equations (l)-( 5 ) to obtain linear growth rate of instability. 
For the inhomogeneous case the equation cannot be Fourier analyzed in x and the eigen 
value can be determined by solving following coupled set of differential equations: 

T a) uPi^nQ^Vi^^y T iP^aXOayi'^OaXlax') 0 (9) 

rji // I / O / I rj-}/ \ f \ I Oq // . 7 Qq, / I 2 I V'' 

rjToaVi^^ + 2-+ TjTq^ Vl^x + -+ rj -+ W 70a - P^anOaVlax 

\ ^Oa / \ ^Oa ^Oa / 

+ioJ-fo^VoayV'i^y = 0 ( 10 ) 

Where rj is ratio of specific heat. For the case of homogeneous plasma, eq. ( 9 ) and eq. (fTOl) 
reduce to 


^IfOaXiay T ^ lOct^lay ^'PaPOa'^lay T 0 (H) 

V'k^OaXiax T ^ 'JOaXlax ^a'^OaXlax T ^^'^OopOayXiay 0 (^2) 
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By taking Fourier transform in x, we obtain the following standard dispersion relation 


(- kl - (^a;^70670p - ^^^7o67op W ^ - r)Tobkl{uj‘^^op - nop) 

V a ^Oa/ V a 


-kl OJ^-iQblQp 


^Oa^Oo 

70a 


nopTiQb y] vl^ + 2nopnofei’opWofe - vTobklriopvlp =0 (13) 


This equation contains two oscillatory mode and one purely growing electromagnetic mode 
which is known as filamentation or Weibel instability. 


III. ANALYTICAL STUDIES 

For analytical tractability of the inhomogeneous problem a choice of sinusoidal variation 
riding on a homogeneous background density such as 

[1 + ecos{ksx)] (14) 

is chosen. Here £ is the inhomogeneity amplitude and kg = 27 rm/L^ (m is an integer) is the 
inhomogeneity wave number with as the system length. To satisfy the quasi-neutrality 
condition the equilibrium beam and plasma density is chosen as 

nob = /d(l -I- ecos(ksx)) 

nop ^ {I -/3){1 + ecos{ksx)) (15) 

where /5 is a fraction. Now choosing the perturbed fields as /„ = 

and assuming e to be small so that retaining only the first order terms (as done in the paper 
by Mishra et al. j^), we can evaluate the growth rate. We evaluate the growth rate Tg,. as 
a f un ction of k, using the same method. The plots shown in Fig. |2] compare the growth rates 
for the homogeneous case with inhomogeneous cases. Fig. |2] (a) , we plot the growth rate F^^ 
versus wave vector k for a homogeneous (solid line) and inhomogeneous cold beam plasma 
system for £=0.1, for kg = 7r(-), 27r(-.-) and 37r(*). The other parameters are nob/nop = 1/9 
or nofe/uoe = 0.1, Vob = 0.9c, Vop = —0.1c, Tob = 0, Top = 0. From this plot we can see 
that for the homogeneous case, the F^^ increases with k for small values of k and saturates 
at /c ~ 3. However, for the inhomogeneous case (£=0.1) Tg^ for small values of k is large 
compared to the homogeneous one. The increase in F^^ with k is very mild and ultimately 
there is a saturation at higher k values. The maximum value of the growth rate F^^ increases 


5 



with kg for a cold system. The effect of transverse beam temperature {TQh± = 10 keV) over 
the growth rate can be seen in Fig. I21[b). The growth rate T^^ of the homogeneous hot beam 
and cold background plasma system (solid line) increases with k for small values of k but 
starts decreasing after k > 0.5 and completely stabilizes at higher wave numbers {k > 1 ). 
The effect of inhomogeneity on the hot beam and cold background plasma can be seen in 
Fig. [21(b) and Fig. [21(c). The Fig. [21(b) illustrate the effect of inhomogeneity scale length kg 
over the F^^ at fixed value of e. We see that for £ = 0.1 and kg = % (which is Xg >c/c(jp), 
the instability domain of F^^ is much larger compared to the homogeneous case. It is also 
observed that when the inhomogeneity scale length is further reduced to half, i.e. kg = 2w{- 
0 -) for £ = 0 . 1 , the instability is suppressed for all values of k. In Fig. [21(c), we plot the 
growth rate Tgr versus wave vector kc/ujp to see the effect of e on Tg^ at Tofej, = 10 keV. We 
see that F^^ reduces with increasing £. Thus the study shows that at finite transverse beam 
temperature the effect of increasing kg as well as £ is to stabilize the Weibel instability. In 
Table - I we have tabulated the value of maximum growth rate for various cases. 

TABLE I 

The maximum growth rate of filamentation instability evaluated analytically under the 
approximation of weak inhomogeneity amplitude. 


%b^{keV) 

£: 

kg 

rp ,.( max .) 

0.0 

0.0 

0.0 

0.2006 

0.0 

0.1 

TT 

0.2057 

0.0 

0.1 

27r 

0.2074 

0.0 

0.1 

Stt 

0.2085 

10.0 

0.0 

0.0 

0.0840 

10.0 

0.1 

TT 

0.0831 

10.0 

0.1 

271 

0.0770 

10.0 

0.1 

Stt 

0.0094 

10.0 

0.2 

Stt 

0.0000 


It is thus noted that for cold beam there is no significant difference between the growth 
rates of homogeneous and inhomogeneous cases except at long wavelengths. However, when 
the beam temperature is chosen as finite the growth rate for the inhomogeneous case sig¬ 
nificantly reduces with increasing kg as well as e. We feel that the reduction in growth rate 


6 






is responsible for transport over long distances observed in targets that had carbon nan- 
otnbes attached to them in experiments. The analytical inference is fnrther corroborated by 
Particle - In- Cell (PIC) simulations. The results of PIC studies are presented in the next 
section. 


IV. PIC SIMULATION DETAILS 


We have performed Particle-In-Cell simulation to study the propagation of beam in both 
homogeneous and inhomogeneous plasma systems. The PIC code PICPSI3D used for this 
simulation has been developed indigenously b y o ne of the authors (Kartik patel) (2]| and 


has been used in several studies in the past 


221 1 . The PICPSI3D was generalized by us 


to study the propagation of electron beam in an inhomogeneous plasma system. Only one 
dimensional variation in space (perpendicular to the beam propagation) has been chosen for 
our numerical studies here, for the purpose of comparison with recently carried out analytical 
linear studies. 


The ions are not allowed to move in the simulations and merely provide a neutralizing 
background for the plasma. This makes the simulation faster and is a valid approximation 
to study the fast electron time scale phenomena of interest. The uniform plasma density Hq 
is chosen to be lOric where ric = 1.1 x 10^^cm“^ is the critical density for Ifim wavelength of 
laser light. The spatial simulation box length = 60 c/cuq, where c/uq = dg = 5.0 x 10“^/rm 
is the skin depth corresponding to the density Hq. The one dimensional simulation box is 
divided in 6000 cells. The grid size is therefore equal to 0.01 dg. Thus, scales shorter than 
the skin depth can be resolved. The total number of electrons and ions chosen for the 
simulations are 1800000 each. This number represents the sum of background nop and beam 
nob electrons. The choice of inhomogeneous ion density and the separation between the two 
electron species of beam and background are made as per Eq. flHl) and Eq. flT5l) 

The beam temperature Tbo± is chosen to be finite in perpendicular direction according to 
Eq.dS]). The time step is decided by the Courant condition. The charge neutrality as well 
as the null value of total current density is ensured initially. The considered system is also 
field free initially as required by equilibrium configuration, the system has an equilibrium 
configuration initially. 
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V. SIMULATION OBSERVATIONS 


For the homogeneous plasma density case (e.g. e = 0.0, kg = 0.0) the system is plagued 
by the usual Weibel instability. This causes spatial separation between the forward and 
reverse electron currents. The separation leads to finite current density in space resulting 
in the growth of magnetic field energy. The evolution of box averaged magnetic field energy 
normalized by Eg = (meCcuo/e)^ energy of system is shown in Fig. 131(b), Fig. 131(b) and Fig. El 
(b). After an initial transient the curve settles down to a linear regime and subsequently 
shows saturation. 

The slope of the linear portion of the main curve has been employed for the evaluation 
of the growth rate of the maximally unstable mode in the simulation. The growth rate has 
been tabulated in Table - II for various cases of parameters. 

TABLE II: 

The maximum growth rate of filamentation instability calculated from PIC simulation. 


Tob±{keV) 

£: 

kg 

Fgr(max.) 

0.0 

0.0 

0.0 

0.2000 

0.0 

0.1 

271 

0.2000 

0.0 

0.1 

Stt 

0.2000 

10.0 

0.0 

0.0 

0.0508 

10.0 

0.1 

271 

0.0349 

10.0 

0.1 

Stt 

0.0333 

10.0 

0.2 

271 

0.0173 

10.0 

0.2 

Stt 

0.0109 


From the table as well as by following the complete evolution it is evident that when the 
transverse temperature of the beam electrons is zero the inhomogeneous density causes no 
change. With increasing value of e and kg the growth rate decreases when the temperature of 
the beam and background electrons is finite. This trend is similar to the behavior of growth 
rate evaluated analytically, shown in Table - I. When the plasma density is homogeneous 
the value of the growth rates evaluated analytically and through simulations are in good 
agreement. However, in the presence of inhomogeneity there is a small disagreement between 






the quantitative values. This can be attributed to the approximate nature of the analytical 
treatment, wherein the inhomogeneity amplitude was assumed to be weak. 

From Figures [3] (a), 0] (a) and [5] (a) it can also be seen that along with the growth of 
magnetic field energy, electrostatic field energy also grows. The development of an electric 
field directed along x during the course of simulation is responsible for this electrostatic 
energy. This electrostatic field develops as a result of the redistribution and bunching of 
electron charges in physical x space. It can be seen from the phase space plots of Figs. Eland 
El that the electrons do reorganize in physical space. Furthermore, the locations where these 
electrons get accumulated are the regions with maximal currents and negligible magnetic 
and electric field as can be seen from Fig. d 

Finally we provide a comparison between the cases of homogeneous and inhomogeneous 
plasmas. It should be noted that the typical scale length of the magnetic field developed 
during the initial phase in the homogeneous case (Fig.lHl(a)) is of the order of the background 
plasma skin depth (e.g. 5.33 x fim). For the inhomogeneous case, the scale length of the 
magnetic field matches initially with the inhomogeneity scale length defined by the choice of 
kg (Fig. El (b)) (provided the scale length of inhomogeneity is smaller than the skin depth) else 
it is determined by the typical value of skin depth. At later stages (the nonlinear phase of 
the instability), however, the magnetic structures coalesce and acquire long scales typically 
comparable to simulation box size in both homogeneous (Fig. |H1 (c)) and in-homogeneous 
(Fig. El (d)) (provided the growth rate remains finite in this case) cases. 

To summarize the main observations are: (i) the inhomogeneous density causes no sig¬ 
nificant difference in the growth rate when the transverse temperature is chosen to be zero, 
(ii) the growth rate of the Weibel instability in the inhomogeneous density case is reduced 
compared to the homogeneous case when the transverse beam temperature is finite, (iii) 
The momentum is typically quite large for beam electrons compared to the background 
plasma electrons, (iv) in the nonlinear regime the typical profile of electrostatic field created 
due to electron bunching in x is similar to that of the magnetic field. The zeros of both 
the fields coincide with each other in space and it is these very locations where electron 
bunching is observed. 
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VI. INTERPRETATION OF NUMERICAL OBSERVATIONS 


We now provide a simplified understanding of tfie observations made by PIC simulations 
fisted out in previous section. In order to understand tfiese results we consider tfie 1-D 
limit (with only variations along x being permitted) of the two fluid system of beam and 
background electrons described in section 11. In 1-D the momentum equations of the two 
electron species from Eqs.(l) are: 


dPxb 

dt 

dpxp 

dt 

dPyb 

dt 

dPyp 

dt 



VybBz 


d(t) , „ ' 


dt 


- VxbB^ 


dt 


^XpBz 


Hb dx 


(16) 

(17) 

(18) 

(19) 

( 20 ) 


Here Pia iox i = x and i = y corresponds to the x and y component of momentum respectively 
for the beam a = 6 or plasma a — p electrons. Also v^a — Piaija (with 7 q. being the 
relativistic factor) is the corresponding velocity. Here Pb represents the transverse pressure 
which is zero for the case when the system is cold. The scalar and vector potentials are 
represented by 4> and A respectively. In 1-D only Ay component is finite. Thus the only 
finite component of magnetic field is along 2 and Bz = dAy/dx. 

The continuity equation can be written as 


duc, 


dt 


^ dric, , 

Wa o A o ® 

dx dx 


( 21 ) 


The Maxwell’s equation become 


dx"^ 


(^STlb 


d'^A, dB^ 


- [ribVyb + npU 


yp\ 


( 22 ) 

(23) 


dB dx 

Here Ub and Up are the total densities of the beam and plasma electrons and 5nb and 5np 
is the difference between the total and equilibrium densities respectively. If one considers 
the transverse temperature to be zero, the linearization of the above set of equations has no 
other term dependent on the electron densities except for Ex = —d4>/dx in the momentum 
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equation. However, Weibel being primarily an electromagnetic instability the electrostatic 
field is very weak. Thus, the predominant term in the momentnm eqnation is dne to sec¬ 
ond term of Vyo^Bz, which is not influenced by the electron density. Thus, in the limit of 
zero temperature the homogeneous and inhomogeneous cases do not show any significant 
difference. 

When the transverse temperature is finite our simulations show reduction in the Weibel 
growth rate. In this case the pressure term in Eq.(13) is effective and depends on the density 
inhomogeneity. It has been shown by an approximate analytical studies in Q repeated and 
presented by ns in Fig.2 that the growth rate indeed decreases in the inhomogeneous case. 
This is the main result of our simulations which qualitatively verifies the approximate results 


of 


0 . 


This can also be relevant for the observed propagation of electrons over long distances 
in the presence of nanowires/nanotubes in experiments [l9(]. 

As we have stated earlier along with the development of magnetic field an electrostatic 
field also develops. This happens due to the bunching of electron densities at the location of 
zero magnetic field as shown in Fig. [3 At the location of zero magnetic field the perturbed 
density shows a maxima and the electric field also passes through zero. This arrangement 
is self consistent. The Lorentz force at these locations vanishes and hence a particle has 
a greater probability to accnmnlate over there. The location of maximnm accumulation of 
electron density in turn results in the vanishing of the second derivative of electric field and 
for a Fourier spectrum this location should correspond to the zero of electric field. 


VII. SUMMARY AND CONCLUSION 

We have shown through 1D3V PIC simulations that the growth rate of Weibel instability 
gets reduced in the presence of density inhomogeneity. This has relevance to a recent ex¬ 
perimental observation of efficient transport of Mega Ampere of electron currents through 
aligned carbon nanotube arrays. The ionization of the carbon nanotubes by the front of laser 
pnlse produces the plasma which has inhomogeneous density. Since the Weibel instability 
gets suppressed in such a inhomogeneous plasma, the current separation is reduced leading 
to the propagation of beam electrons over large distances. 

This mechanism of efficient electron transport was earlier invoked by Mishra et al. 0 
wherein it was shown analytically nsing two flnid description that the Weibel instability 
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gets suppressed. The present work supplements it with PIC studies. Our PIC simulations 
support the analytical observations qualitatively. The quantitative values of the growth rate 
differ slightly showing the approximate nature of the analysis. 

Acknowledgement: We thank S. Mishra for many useful discussions. 
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FIG. 1: schematic of model configuration 



FIG. 2: Growth rate of the ID filamentation instability versus the wave number k for a 
homogeneous and inhomogeneous beam-plasma: (a) shows growth rate for cold 
homogeneous beam-plasma system (solid line), for parameters Uofe/noe = 0.1, VQt, = 0.9c, 
Vop = —0.1c and for inhomogeneous, £=0.1 at ks=7r(-), ks=7r (-.)and ks=7r{*). (b) at 
transverse beam temperature T;,o±=10 keV, homogeneous (solid black), inhomogeneous for 
£=0.1 at ks=TT{- -), 27r(-.) and 37r(-*) (c) at transverse beam temperature Tbo±=10 keV, 
homogeneous (solid black), inhomogeneous at ks=37r for £=0.1 (- -) and 0.2 (-.) 
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FIG. 3: Temporal evolution of the field energy densities for cold homogeneons and 
inhomogeneous beam-plasma simulation, (a) normalized electrostatic x-component of 
electric field energy (b) normalized z-component of magnetic field energy 




FIG. 4: Temporal evolntion of the field energy densities for hot (T^ox^lO keV ) 
homogeneous and inhomogeneous beam-plasma simulation, (a) normalized electrostatic 
x-component of electric field energy for £=0, /Cs=0 and £=0.1 and 0.2 for ks=27r (b) 
normalized z-component of magnetic field energy for £=0, ks=0 and £=0.1 and 0.2 for 

ks=2Tr 
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FIG. 5: Temporal evolution of the field energy densities for hot (rfco±=10 keV ) 
homogeneous and inhomogeneous beam-plasma simulation, (a) normalized electrostatic 
x-component of electric field energy for e=0, ks=0 and £=0.1 and 0.2 for ks=3Tr (b) 
normalized z-component of magnetic field energy for £=0, ks=0 and £=0.1 and 0.2 for 

ks=3TT 
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FIG. 6: projection of f (x ,Px) for homogeneous cold beam plasma at u)ot= 25.1789 and 

74.0555 
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FIG. 7: Bunching of perturbed electron 
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FIG. 8: spatial configuration of normalized magnetic field homogeneous and 
inhomogeneous hot beam plasma (a) £=0.0, ks=0 at a;ot=49.3703 (b) £=0.2, ks=37: at 
070^=49.3703 (c) £=0.0, ks=0 at C(Jot=74.0555 (d) £=0.2, ks=3'K at o;ot=74.0555 
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FIG. 9: projection of f (x ,Px) for warm system (a) beam with £=0.0, fcs=0 at C(Jot=49.3703 

(b) beam with £=0.2, ks=3Tr at a;ot=49.3703 
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